Glass transition of hard spheres in high dimensions 
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We have investigated analytically and numerically the liquid-glass transition of hard spheres 
for dimensions d — >■ oo in the framework of mode-coupling theory. The numerical results for the 
critical collective and self nonergodicity parameters fc {k; d) and fi"^ (k; d) exhibit non-Gaussian k- 
dependence even up to d = 800. /i"' (fc; d) and fc {k; d) differ for k ~ d^^^, but become identical on 
a scale A; ~ d, which is proven analytically. The critical packing fraction (pc (d) ~ d?2~'^ is above the 
corresponding Kauzmann packing fraction tpK (d) derived by a small cage expansion. Its quadratic 
pre-exponential factor is different from the linear one found earlier. The numerical values for the 
exponent parameter and therefore the critical exponents a and b depend on d, even for the largest 
values of d. 

PACS numbers: 64.70.P, 64.70.Q, 64.70.kj 
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I. INTRODUCTION 

In many situations the analytical treatment of a spe- 
cific physical problem simplifies drastically if the spatial 
dimension d becomes infinite. For instance, it is well- 
known that the mean field theory for systems in thermal 
equilibrium becomes exact for d = oo. The equation of 
state of a fluid can be obtained from a virial expansion. 
For the fluid of hard spheres it has been shown that for 
such packing fractions ip for which the second virial term 
(which is proportional to tp) is finite or at most alge- 
braically increasing with d, the third and higher order 
virij,l terms vanish exponentially fast in the limit d ^ oo 

Hard sphere systems are ideal systems to study not 
only equilibrium properties, but also the liquid-glass 
transition and glassy dynamics. The mode-coupling the- 
ory (MCT) Q is a microscopic theory of an ideal glass 
transition. Knowledge of the static density correlators al- 
lows to calculate the long time relaxation of a supercooled 
or supercompressed fluid and to locate the glass transi- 
tion point at which the ergodic behavior in the fluid phase 
changes discontinuously into a nonergodic one. The ex- 
perimental results for colloidal fluids, which can be mod- 
elled by hard spheres, exhibit agreement with the cor- 
responding MCT result after the transient regime over 
several decades in time within ten percent precision 0, Q . 

Properties of equilibrium phase transitions, e.g. the 
critical exponents at a second order phase transition, de- 
pend strongly on the spatial dimensionality. This has mo- 
tivated the investigation of the glass transition for d = 2 
P-Q and d = 3,4 [8l-[To|. The most important approxi- 
mation of MCT is the factorization of the memory kernel 
Q. This kernel is a time- dependent four-point correlator 
of the density modes p{k) which is approximated by a 
product of time dependent two-point correlators. This 
factorization resembles the mean field approximation re- 
placing a static two-point correlator, e.g. the spin-spin 
correlator (SiSj) for an Ising model, by a product of the 
order parameter, e.g. the magnetization (Si) in case of 
the Ising model. Based on this analogy, MCT has been 



interpreted as a mean field theory with the two-point den- 
sity correlator as an order parameter [111] , where spatial 
fluctuations of the correlation between the pair densities 
p{f,t)p[r-\- S,t -\- t) and p{f',t')p{f' -\- S,t' -\- r) are ne- 
glected. In a next step, these spatial fluctuations are 
taken into account. Finally it is shown that the upper 
critical dimension where the spatial fluctuations do not 
influence the critical behaviour is c?c = 6 [ll|, for sys- 
tems without and dc = 8 [isj with conserved quanti- 
ties. This implies that the square root singularity of the 
nonergodicity parameter and the relation between the 
exponent parameter X(d) and the "critical" exponents 
a(d), b{d) i are universal above dc However, the 

exponent parameter X{d) itself being determined by the 
static structure factor at the glass transition singularity, 
depends on d. The interpretation of MCT as a mean field 
model challenges the investigation of MCT with full k- 
dependence for d — oo. As already mentioned above, an- 
alytical calculations simplify for d — > oo, e.g. the leading 
order term of the static and direct correlation function 
for hard spheres are known and become dominant (see 
below). Consequently we will focus on the MCT glass 
transition of hard spheres in high dimensions. 

Let us shortly review what is already known for hard 
spheres and d — >■ oo. Taking for the direct correlation 
function the leading order of a virial expansion (see be- 
low), using the Vineyard approximation [l5| for the nor- 
malized collective nonergodicity parameters f{k;d), i.e. 
it is /(fc; d) « (^; ^-iid assuming the nonergodicity 
parameters f'^^\k\d) of the self correlator to be Gaus- 
sian in fc with width a, a self consistency equation for 
a follows from MCT. As critical packing fraction for the 
glass transition it has been found 16] 
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The replica theory for the structural glass transition 
p7j is another microscopic theory. It allows to calcu- 
late the Kauzmann temperature Tk or the corresponding 
packing fraction LpK at which the configurational entropy 
per particle vanishes. Applied to hard spheres in high 
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dimensions and performing a small cage expansion it is 
found pJi] 



coordinates reads: 



ipK {d) = din (d/2)2~ 



(2) 



Our main motivation is to explore the MCT scenario 
for d — >■ oo, i.e. we want to investigate whether the A2 
singularity Q of MCT in d = 3 survives for d — ?► 00. 
Furthermore, we want to check whether the critical non- 
ergodicity parameters /c(fc; d) and /i'*'' (fc; d) of the collec- 
tive and self correlators, respectively, are Gaussian and 
become equal for d — )■ cxd and whether the critical pack- 
ing fraction Lpc{d) coincides with (pf-^ (d) and if not how 
it relates to the Kauzmann value ipxid). In a first step 
we have solved numerically the MCT equations for the 
collective and self nonergodicity parameters /(fc; d) and 
/('^^(fc; d), respectively up to d = 800. Inspired by the nu- 
merical solution we have investigated in a second step the 
corresponding equations analytically. The outline of our 
paper is as follows. The next section presents the MCT 
equations in arbitrary dimensions d and their numerical 
solution for the nonergodicity parameters of the collec- 
tive and self correlator. Based on these numerical results 
we present in the third section an analytic investigation 
of the MCT equations for hard spheres for d — >■ 00. The 
final section IV contains a summary and conclusions. 



II. MCT EQUATIONS AND NUMERICAL 
SOLUTION 



A. MCT equations 

We consider N hyperspheres with diameter tr in a d- 
dimensional box with volume V. The number density is 
n = N/V and the packing fraction: 



if = nVd (a/2) 



(3) 



k+p 



dp / dqV {k,p, q) <j) {p,t) <j) {q,t) 

J\k-p\ 



with the vertices in arbitrary dimensions d 
V {k,p, q)^n -j^S [k) S {p) S {q) ■ 



(6) 



4fc V - (fc' + - q^) 



((J-3)/2 



(7) 



■[{e+p^-q')c{p) + {k'-p'+e)c{q)\ 



c{k) is the direct correlation function and fid = 
27r'^/^/r(d/2) is the surface area of a d-dimensional unit 
sphere. 

The corresponding equation of motion for the self cor- 
relator i/i'^^ (fc, t) follows from Eq. ([5]) by replacing 7^ and 
m{k,t) by 7^'*-' and rrS^^kjt)^ respectively. m^'''\k,t) is 
given by: 



m(^) {k,t)^4''^ iq,t) 



k+p 



(47r)" 



dpi dqV^'Hk,p,q)cj,{p,t)cj,^'HQ,t) 
J\k-p\ 



(8) 



with the corresponding vertices Q: 



(fc,p, q) = 2n^S (p) [{k' +p'-q^)c( 



Aep'-{e+p'^q^Y 



(9) 



Note, that the vertices Eq. JT]) and Eq. ^ reduce to 
the well-known expressions [3| for d — 3 for which the 
triple direct correlation function c^'^^{k,p, q) has been ne- 
glected. 
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B. Static correlation functions 



Vd (R) = 



rd/2 



r(d/2-t-i) 



(4) 



the volume of a d-dimensional sphere with radius R. T{x) 
is the Gamma function. 

MCT provides an equation of motion for the intermedi- 
ate scattering function S{k,t) Q. For a one-component 
liquid with Brownian dynamics the MCT equation for 
the normalized correlator cj){k,t) — S{k,t)/ S{k) is given 
by: 



-/k(j)ik,t) + (j){k,t)+ / dt'm{k,t-t')(l){k,t') ^Q. (5) 

JO 

7fe is a microscopic relaxation rate related to the short 
time diffusion constant. The memory kernel in bi-polar 



The static correlation function S{k) = S{k,t = 0) is 
related to the direct correlation function by the Ornstein- 
Zernike equation: 



S {k; d, (p) = [I — n (ip) c (fc; d, (p)] 



(10) 



The direct correlation function c(fc; d, ip) is known ana- 
lytically for d — > 00, in case that the third and higher 
virial terms of the virial expansion can be neglected. It 
has been shown [l^ that the truncation for d — > cx) at 
the second virial term is even valid above the packing 
fraction at which the virial series diverges. Under these 
conditions it is c{r;d,(p) = —Q{(j — r) — f{r) (Mayer 
function) from which one obtains 

c{k; d,^)-c (fc; d) = - {27rf^ a''Jd/2 (ka) / {kaf^ 

(11) 
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where c{k; d) does not depend on (p. a is the diameter of 
the hard spheres and J„(x) the Bessel function of order 
n. Note, that the d- and (^-dependence of the various 
quantities is made exphcit in cases where it is useful, 
and suppressed otherwise. There are two d-dependent k- 
scales on which the fc-variation of S{k; d, ip) is different. 
For 



k = ka/\/d, (p ^ 



(12) 



it follows from Eqs. ([3]), ()10|) and by using the Taylor 
series for Jd/2(\/rf k) [20] in the scaling limit k — >■ oo, 
d — >■ oo, (/?—>■ such that k — ka j \fd and (p — 1'^ip are 
fixed: 



jmi ^((yd/o-j fc;d,2-''(^j = \V-pc{k)\ ' = S {k;<p) 

(13) 

where: 



exp ( -^fc^ 



(14) 



Figure [T] demonstrates the convergence of S (fc; d, p) to 
the scaling function S{k\ p) on the scale k for (p — 2'^p 
fixed. 
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FIG. 1: (colour online) S (^(^/cr)fc; d, 2-''(^) on the scale k 
for = 1 and d = 5, 10, 20. The bold black line is S (k; l). 

S(k; (p) does not exhibit peaks. It increases monotoni- 
cally from S{0; p) = 1/[1 + 1^] < 1 to the ideal gas behav- 
ior S{k] <p) = 1 for k — > 00. This is the well-known effect 
for /c — > of the suppression of the compressibility below 
the corresponding value of an ideal gas which, however, 
is much weaker than for d = 3. 

The second /c-scale is linear in d: 



fccr/d, <p = p2'^e-'^l'^ 



(15) 



Making use of the asymptotic expansion of Jn{nx) [2C 
one obtains: 



S ({d/a)~k;d,2-'^e'^/^ip 



1 - ipcd k 



Sd[k;p 
(16) 



where 



-1/4 



• exp 



arctanhA/ 1 — 4fc 



(17) 



l-4/c2 



for fc < 1/2, and 

(i). 



Cd 



■ cos 



Ak^ - 1 



-1/4 



(18) 



- J 4fc2 - 1 - - arctan J Ak"^ - 1 - - 
2 V 2 V 4 



for k > 1/2. Sdik-ip) oscillates for k > 1/2. 

The position fc* (d) of the main peak (first sharp diffrac- 
tion peak) of S{k; d, ip) is given by the first nonvanishing 
zero of Jd/2+i{^)j which is [20| : 

fc, (d) = (d/2 + 1) + ao (d/2 + 1)^^^ , oq = 1.8557571. 

(19) 

Since /c* (d) and c(fc; d) from Eq. (fTTI) are (^-independent 
the packing fraction (/3*(d) for which S{k^,{d);d,(p) di- 
verges follows from Eqs. (g]) and (fTOl) : 



(d) 



^''/^ ((7/2)" 



F(d/2 + l)c(fc* (d);d)' 



(20) 



Using again the asymptotic properties of the Gamma and 
Bessel function and especially from 



Jd/2{K{d)) = J'd/2+iihm = -bo-{d/2 + n 
bo 1.1131028 



-2/3 



(21) 



(d) 
Co 



Co • d^/^ exp (ao (d/2)^/^) (^8^ 



&^i7r-i/22-2/3 



0.867956, 



(22) 



i.e. the leading d-dependence of (/?* (d) is the expo- 
nential factor (v/s/e)-'* (1.7155)"'^ [2l|. Note that 
S{k;d,p>) > for all k provided (p < (p*{d). For 
p> <^ p>* (d) the height of the first sharp diffraction peak 
of S{k;d) is given as 



Sik,{d);d,^) ^ 1 + 



P* (d) 



(23) 



Figure [2] presents S{k;d,p)) on the scale k = ka/d for 
d = 200 and ip of order 1^9* (d). 

S{{d/a)k; d, (p) is practically zero for k < 0.5, develops 
its main peak at fc*(d) ~ k^{d)a/d and decays very fast 
to one for k > fc*(d). In contrast to this. Figure [3] shows 
S{{d/a)k;d,p>) again on the scale k but for p> of order 
ipc{d) = 0.22 •d^2~'^, which will turn out to be the critical 
packing fraction of the MCT glass transition. Note that 
p>c (d) is exponentially smaller than pi^ (d) . 
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FIG. 2: (colour online) S {{d/a) k; d, if^ versus k for d — 200 
and ip — O.lift, (p = O.Sif,. 
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FIG. 3: (colour online) S ( {d/a)k\ d, ip ) versus k for d = 200, 
400, 800 and ip ^ ifc {d). 



Except for k — 0{l/Vd) which is ka = 0{^fd) it is 
S{k]d,Lp) = 1, i.e. the static structure factor on a k- 
scale linear in d is very close to that of an ideal gas for 
d — 00. The first sharp diffraction peak of the conven- 
tional liquids has disappeared due to the lack of interme- 
diate range order at ip = (pc{d). These results will be used 
for the analytical treatment of the MCT equation. It is 
important to note that [S{k; d, ip) — 1] is very small for 
ka = 0{d), (fi = 0{fc{d)) and d ^ 1, but not zero. Ac- 
cordingly the direct correlation function does not vanish, 
in contrast to an ideal gas with finite density. Therefore 
the vertices Eq. ([7]) and Eq. ([9]) are nonzero and exhibit 
nontrivial fc-dependence. 



C. Numerical solution 

The nonergodicity parameter for the collective correla- 
tor is the long time limit of the normalized intermediate 



lim (p {k, t; d, ip) 



scattering function 

/ {k;d, ip) 

and similarly for the self correlator 

/(^)(fc;d,^)= lim ^('^J {k,t;d,^). 



(24) 



(25) 



They are the order parameters for the liquid-glass tran- 
sition. From Eq. ([S]) and the corresponding equation for 
(ft^'^^kjt; d, (p) one obtains the algebraic, nonlinear equa- 
tions for the nonergodicity parameters 

/ (fc; d,ip)/[l-f (fc; d, ip)] = Tk [f (q; d, ^) ; d, p,] (26) 

and a similar equation for f^^\k;d,ip) by replacing J-^ 
by J'j, ' . Note, that we also made explicit the d- and tp- 
dependence of the functional J-k on the r.h.s. of Eq. (|26l) . 
Eq. ([25)) and the corresponding one for f'^'^^ {k;d,p) has 
been solved numerically as follows. 

Eq. (1251) is rewritten such that the nonergodicity pa- 
rameters / (fc; d, pi) can be evaluated by iterating the 
equation 



(fc;d,^) = 



Fk [/W {q;d,p)-d,p>]+l 
with the initial value 

/(")(fc;d,^) = l 



(27) 



(28) 



and similar equations for /^*-' (fc; d, (p). Note, in case 
of hard spheres the functional J-^ for the zero or- 
der iterate /*^°-' (fc; d, (/?) from Eq. ([25]) exists only for 
a finite cut-off at fcmax- The integrals appearing in 
[Z*-*^ (fc; d, p) ; d, p\ are replaced by Riemann sums 
with an upper cutoff crfcmax = max(40d-^/^; 4d; 0.2d^/^) 
and 500 gridpoints for d < 200, 1000 gridpoints for 
200 < d < 600 and 1500 gridpoints for d > 600. The crit- 
ical packing fraction Pc{d) is the packing fraction, where 



/(fc;d, p) 



= , (y9 < pc{d) 
, p> Pc{d) 



(29) 



and the critical nonergodicity parameters are given by 

f,{k;d) = f{k;d,p,id)). (30) 

Because the real critical packing fraction and the critical 
nonergodicity parameters can never be computed numer- 
ically in finite time, we evaluated /c (fc; d) at a packing 
fraction ipdd) where limi_j.oo Z*-*^ (fc; d, ipc{d)) = but 



min ( max 

i \ fc 



(fc;d,(^e(d))-/« (fc;d,^,(d)) 



{k;d,p,(d)) 



< e 



(31) 



with e = lO"'^ for d < 600 and e = 10"^ for d > 600. 
It can be estimated that the relative difference between 
this packing fraction ipdd) and the real critical packing 
fraction pc{d) is of order e. It has been verified that the 
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system really becomes nonergodic near this packing frac- 
tion, i.e. / (fc; d, (/?) ^ for > (1 + e • 0{l))^c{d). The 
critical nonergodicity parameters can then be approxi- 
mated by fc{k;d) = f^^°^ {k;d,ipc{d)) where iq equals 
the iteration step, where 



max 

k 



/(■'o+i)(fc;d,(^,(d))-/(^o)(fc;d,(^,(d)) 



(32) 



reaches a minimum |24| . It has been verified that there 
are no visible differences in the critical nonergodicity 
parameters obtained by this procedure with different 
values of e and that f{k;d, ipc {d) + Aip) converges to 
/(*°-' {k; d, ipc{d)) with order y^Aip. Additionally it has 
been verified that fc (fcinax;^) < 10""'^^ for all evaluated 
dimensions. By increasing ftmax and the number of grid- 
points the relative error of the critical packing fraction 
due to the discretization can be estimated to be smaller 
than 10-3 for d < 600. 

The nonergodicity parameters always show numerical 
artifacts on the first few gridpoints in fc-space. This is 
a problem when trying to observe the characteristics of 
fc{k;d) for small wavenumbers, especially for high di- 
mensions. So we interpolated fc{k;d) onto a much finer 
fc-grid and performed one single iteration step equivalent 
to the one given in Eq. (I27|) . This procedure improves 
the result for fc (fc; d) by shifting the numerical artifacts 
to much smaller values of fc. 

From this solution we obtain the critical packing frac- 
tion ipc{d), shown in Figure 



10^ 



■D 



10^ 



10 



10' 



10^^ 
10 



100 

d 



1000 



FIG. 4: (colour online) d-dependence of the critical packing 
fraction kpc (d) on a log-log representation. The dashed line is 
2''ipc{d) ^ ad^. The dotted line is (ys* (d) from Eq. pj)) . 

The d-dependence of (pc can be well fitted by 

^c{d)^ad^2-'^, 0^^ 0.22. (33) 

The critical nonergodicity parameters /c(fc; d) and 

Jc^\k\ d) are presented in Figure[S^ and [SJd, respectively, 
for different values of d. 

Figure [S^ and [S}d reveal the following properties: 
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FIG. 5: (colour online) fc-dependence of the critical nonergod- 
icity parameters fc (fc; d) (a) and (fc; d) (b) for d — 200, 
400, 600 and 800. Diamonds in (a) mark the numerical pre- 
cise values for fc (0; d) and the dotted line in (b) presents 
/c(fc;600). 



(i) fc{k;d) and fc'^\k;d) differ on the scale ka = 
0{Vd), but become identical on the scale ka = 
0{d), for d large enough. 



(u) Both, /c(fc;d) and /, 
fc-dependence. 



{k;d), exhibit non-Gaussian 



(iii) There are three characteristic fc-scales. fc{k;d) in- 
creases from fc{0;d) to its maximum value on a 
scale ka ~ Vd, develops a plateau on a scale ka ^ d 
and it varies on a scale ka ~ c?'^/^, on which a 
steep descent to zero occurs for ka ^ kod^^'^ where 
fcg = 0.15. The plateau value on scale ka d con- 
verges to one, for d — c«. 

(iv) Since fc {k;d) changes from one to zero around 
ka = fcod'^/^ we will choose fcg such that 
fc(kod^'^/a\d) = 1/2. fc{k;d) is of order one for 
{ka - of order d^'^, i.e. for (fc - kod^'^) of 
order c?~^/^. 

Using Eq. (|26l) . fc{0;d) can be represented as a func- 
tional of fc{k;d) Q]. Making use of this relationship 
yields the numerical precise values for /c(0; d) shown by 
diamonds in Figure [5^. Note that for the self correlators 
it is /i*' (0; d) = 1 for all d, because the momentmn of a 
tagged particle is not conserved. 

A crucial quantity of MCT is the exponent parame- 
ter A(d) which determines the critical exponents of both 
power laws in time, the critical law and the von Schwei- 
dler law, and the divergence of the corresponding relax- 
ation time scales at the glass transition singularity [sl]. 
Figure [5] presents the numerical result for A(d) with an 
estimated relative error of about 5 • lO"'^ for d < 600 and 
possibly a larger error for d = 700 and 800. 
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FIG. 6: d-dependence of the exponent parameter A. The full 
symbols mark the regime where (pc{d) from Figure [4] follows 
the asymptotic result (|33p . The values for d = 700 and 800 
depicted by full circles possibly have a larger relative error. 



Since the direct correlation function pT|) is not correct 
for small d, the variation of A with d below 100 and par- 
ticularly the high sensitivity close to d = 17 is an artifact 
of the incorrect static input. This holds also for <Pc{d) of 
Figure SI The concave curvature of ipc (d) and the cusp- 
like behavior of A (d) around d — 17 relates to the fact 
that the glass transition is still influenced by the main 
diffraction peak of the static structure factor for d < 17, 
while this peak is not important any more for d > 17. 
Figure H] reveals the correct asymptotic d-dependence to 
appear for d > 100. 



III. MCT EQUATIONS: ANALYTICAL 
APPROACH 



In this section we will demonstrate that the MCT func- 
tional J-k strongly simplifies for d — )■ oo. The essential 
steps will be given, only. 

In a first step we rewrite Tk[f{q)] (Eq. (O with (t){q,t) 
replaced by the nonergodicity parameters f{q)) on the 
scale k = ka/d. Quantities on this scale will be denoted 
by a tilde, e.g. /(fc). Now we prove that F[f{q)] on this 
scale reduces to J^^*'* for k of order one and larger. 

The last square bracket in Eq. ([7]) contains a mixed 
term ^ c{pd/a)c{qd/a) which is oscillating in p and 
q faster and faster under an increase of d. Since 
S{pd/a) S{qd/a) and f {pd / a) f {qd / a) are smooth and 
not strongly oscillating (see section II) the mixed term 
integrated over p and q will not contribute for d — > cx). 
Taking account of this fact and using the integrand's 
symmetry with respect to p O q we get for d ^ oo the 



MCT functional Eq. ^ but with the replacement: 
V (jcd/a,pd/a,qd/aj — > 
^ S {kd/a) S (qd/a) V^''> {kd/a, pd/a, qd/a^ 

V^""^ (kd/<7,pd/a,qd/ay (34) 

The latter step in Eq. uses S{kd; d) 1 for d ^ oo 
and k of order one or larger. Remember that this does 
not hold for k = 0{1/Vd) (see section II. B). Eq. (|34l) . 
together with the Vineyard approximation for (j) (pd/a, t) 

in J^l''^ (Eq. dH])), implies that the MCT equations for 
the self and collective correlator on a fc-scale linear in 
d become identical for d — > oo. However, for large but 
finite d, there is an interval ka G [0, 0{^/d)] for which 
both MCT functionals differ from each other. 

Having reduced V to y'^-* on the fc-scale linear in d 
further simplifications occur due to d — oo. First of 
all we can replace S (pd/cr) in V^'^^kd/a, pd/a, qd/a) by 
one. The product of both square brackets in © can be 
rewritten as 



-2d+2 j2d-2 



2kp 



k"^ +p^ - q^ 

2kp 



iPd/a) 



The square bracket in 
e ^ V 2fcp J ^ Then we use 




X e 



3/2 




which allows to perform the integration over q. Next 
we account for the asymptotic behavior of J„ (nx) for 
n — ^ oo [IQ] and obtain from Eq. (fTTj): 

[pd) = {2T:)%'''d-''j%, [{d/2) 2p] 
4 {2T:f a^'^d^'^ [p'^Tidy/Ap^ - l) 6 (^p - ^ 



— a/ Afp — 1 arctan \/ AfP — 1 

2 2 A 



(37) 



With these simplifications and the fact that cos^ [. . .] in 
(|37| oscillates very fast for d large with average 1/2 we 
arrive at: 



f{q) 



k^nd 



dp 



P 



fip)ifiq-) + fiq+) 



where: 



q± 



p^±2\ 



1/2 



(38) 



(39) 
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Our first goal will be the evaluation of the critical pack- 
ing fraction {d). For this we choose k = ko = kod^^'^ 
such that /c(fco) = 1/2 (see section II. C). Eq. (^5)) im- 
plies 



fc (q) 



= 1, 



if ^ ipc (d) . 



(40) 



fc{Q±) is of order one for k — ko and p = O (1) and 
decays rapidly to zero for p 3> 1 . Furthermore / (p) ^ 1 
for p = 0(1). Therefore the integral in Eq. ([55)1 is of 
order one, i.e. order d'^. Then we get from Eq. (l38l) with 

tp = ipc (d) and k = k^/d: 

1 = -^So [/c (9)] = Vc (d) (2Vd') fco^'O (d°) . (41) 

kf) is of order d^\ as well. Consequently it must be: 

(p^ [d) = const ■ d^2^'^ (42) 

in agreement with the numerical result for c? > 100 (cf. 
Figure S]). Next we choose k < ko. fc{q±) ~ for 
q± > ko, which implies 



'{2/d) F 



k"^ — fc2 



for d ^ oo and k < ko. With fdp) { fdq-) + fc{q+) 
2 for p < 



fcg — fc2 the integration in Eq. (155)) can be 

done. Then we get from psp for ko — k ^ ko/Vd after 
substitution of (rf) from Eq. ([33)) : 



lim ( J- 



(-^v^fc [/c (?; c?) ; d] /Vd) = Fo (fc) (43) 



with the master function 

-1 



^^0 fc 



an " fc ^ A / fcj? — fc2 , fc < fco 



, fc > fco 



(44) 



Figure [7] presents the numerically exact result for 

■^Vdk -^c ^) ' ^ /v^ as function of fc = ka/d^^^. 

This figure demonstrates the convergence of ^ ^/dk/ 
at the glass transition singularity to the master function 
Fo(fc)- For fc < fco, i.e. fc < fco the critical nonergodicity 
parameters /c(fc; d) are close to one for d oo. Making 
use of Eqs. ()26p and (H5)) the fc- and d-dependence of 
/c(fc; d) can be expressed as follows 



, , VdFo (k/Vd) 

/c(fc;rf = ^} . 

^ ^ 1 + VdFo [k/Vdj 



(45) 



i.e. on the scale fc = ka/d^/"^ it is 



lim 



/c ({d'"^/<T)k d) = (fc) = e (fco - fc) . (46) 



The convergence of the critical nonergodicity parameters 
to a step function is demonstrated in Figure |H| 




0.001 



FIG. 7: (colour online) The numerical result for 
^Vdk [/<= (^i c') ',d^ l^fd as function of k for various d. The 
master function Fo{k) is shown as bold black line. 
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FIG. 8: (colour online) The numerical critical nonergodicity 
parameters on the scale k = ka/d^''^ for various dimensions 
d. 



IV. SUMMARY AND CONCLUSIONS 

The liquid-glass transition for hard spheres in high di- 
mensions d has been reinvestigated in the framework of 
MCT. Our aim has not been exploring the validity of the 
MCT approximations for d — oo (we come back to this 
point below) but to take MCT as a microscopic theory in 
any dimension and to check the generic MCT-bifurcation 
scenario {A2 singularity) and to calculate the critical 
packing fraction, the corresponding nonergodicity pa- 
rameters and the exponent parameter. The direct corre- 
lation function for hard (hyper-) spheres for finite d is not 
known exactly. However, if ip — 0{ipc{d)) ^ d^2~'^, being 
exponentially smaller than Lp^,{d) at which the Kirkwood- 
like instability occurs, the corrections to the leading order 
result fEq. [TT)) can be neglected for d — > cxd. This offers 
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the possibility to calculate various quantities like (pc{d), 
fc{k;d), f^'\k;d) and \{d) from MCT for a liquid of 
hard spheres in high dimensions. 

A. Summary 

Let us first summarize our results. The numerical so- 
lution of the MCT equations for the collective and self 
nonergodicity parameters up to d = 800 reveals non- 
Gaussian dependence of the critical nonergodicity pa- 
rameters /c (fc) and /i*'' (fc) on the wavenumber k (cf. 
Figure [5]) . Three different /c-scales were found on which 
fc{k) behaves differently. For ka = 0{Vd), fc{k) in- 
creases from fc (0) < 1 to a maximum value close to 
one. fc (k) stays close to one for ka = O (d) and fi- 
nally drops to zero for ka = k^d^/"^. The decrease to 
zero happens on the scale k = ka/d^/"^ in an interval 
around fco with width of order 1/d . fc^\k) and fc{k) 
are identical for ka = 0{d) but differ for ka = Old}^^). 
The exponent parameter X{d) (cf. Figure |6l) varies with 
d even for d > 100 and strongly exceeds the value 
A(3) =2; 0.735 |25i] for d > 1. Note that \{d) can- 
not be larger larger than one. The critical amplitudes 
h{k] d) [s'] exhibit the expected fc-dependence. They are 
in anti-phase with fc{k;d), i.e. they have a maximum 
(minimum) where fc{k;d) has a minimum (maximum). 
Particularly, on the scale ka — 0{d) the critical ampli- 
tudes are rather small since fc{k]d) ~ 1. The numer- 
ical results (up to d = 800) also show that the largest 
eigenvalue Eo{ip) of the stability matrix [s*] is not degen- 
erate and it approaches the bifurcation point at ipdd) 

as (1 - ^^o('^)) ~ Viv- Vc{d))/ipcid), for ip ipc{d) 
from above. Hence the glass transition singularity is 
an A2 singularity, consistent with an exponent param- 
eter X(d) smaller than one, as demonstrated by Figure 
m What happens ior d — 00 is not clear. The distances 
between the largest eigenvalues of the stability matrix 
obtained numerically slightly decrease with increasing d. 
Whether or not the largest eigenvalue becomes degener- 
ate for d = cx) is an open question, similar to the question 
whether X{d) in Figure [5] converges to one for d 00, 
or not. Therefore, our numerical results do not allow to 
exclude a higher order singularity at d = 00. 

Inspired by these numerical results we have been able 
to prove analytically that the Vineyard approximation 
becomes exact for d — > cx) on a scale ka = 0(d), and 
that the critical packing fraction ipc (d) decays exponen- 
tially as 2^'*, with a prefactor which is quadratic in d. 
This is consistent with the numerical result for d > 100. 
Furthermore, the analytical approach has also shown that 
the fc- and d-dependence of the critical nonergodicity pa- 
rameters and of the MCT functional J^j;.[/c (9; d) ; d] at 

ifc (d) can be obtained from a master function Fo^k) (cf. 
Eqs. (1121) and (gS])). This relationship yields 

^lim fc {ka-^d^'^\d^ = 6 (fco - A;) , d ^ 00 (47) 



where fco — 0.15. 



B. Validity of MCT 

Although it has not been our purpose to prove or dis- 
prove the validity of MCT for d — >■ 00, it might be useful 
to comment on this question. First of all, the vertices ([7]) 
and seem to be exact for d — > 00, since the leading 
order of c(fc; d, ip) is known analytically and the neglec- 
tion of the triplet direct correlation function c^^^{p,q), 
which also enters into the vertices |3|], is justified for 
(f ~ d^2~'^ and d — cx) (see appendix). Accordingly 
the convolution approximation for the static three point 
correlator {p{k)* p {p) p {q)) (usually done in MCT) be- 
comes exact for d ^ 00 and packing fractions such that 
2'^(p{d) does not increase exponentially or faster with d. 
The factorization of the static four-point density corre- 
lator S'(^)(9i,92,g3,94) = jj {p* {qi) p* {q2) piqa) p{q4)), 

which is needed for the projector onto pairs of density 
modes, is another approximation. Similar to S^^'^{p,q) — 
jf{p{p + q)*p{p) p{q)) (see appendix) one can define a 
quadruplet direct correlation function S-^^ (91,92,93) via 
a corresponding Ornstein-Zernike equation. However, 
this equation is already rather involved {261] such that 
we have not attempted to prove that the factorization 
is valid for ip — ip^ (d) and d — >■ cx). So it remains an 
open question whether this factorization becomes exact 
again. If so, the "only" two remaining crucial approx- 
imations of MCT are the projection of the fluctuating 
force onto pair modes and the subsequent factorization 
of the pair density correlator with reduced dynamics into 
a product of density correlators (fc, €) with full dynam- 
ics. Whether these two steps become exact for d — > 00 is 
an interesting but also a highly nontrivial question. Acti- 
vated processes smear out the glass transition singularity. 
Since it seems that the local barriers between adjacent 
metastable configurations increase with increasing d (see 
also Ref. [18]) these hopping processes may become sup- 
pressed for d — >■ cxD. Even if this is true, it is not obvious 
that the remaining two approximation of MCT become 
exact for d — ^ cxd. 

One might conclude that MCT for d 00 necessarily 
fails because of S(k\d,pic{d)) = 1 for ka — 0{d), exclud- 
ing the cage effect [3| as driving mechanism and one may 
argue that the dynamics will be described better by a 
Boltzmann-Enskog equation. Indeed, a modified Enskog 
equation was derived from the binary-collision expansion 
[27| . But its validity for packing fractions of order 2"*^ 
has not been proven. Here, two comments are in order. 
First, a cage does not necessarily require a maximum 
number of adjacent spheres. For instance a (simple-) 
hypercubic lattice built up of periodically arranged hy- 
perspheres has 2d nearest neighbors and a volume frac- 
tion fhypercubeid) ~ d~'^^^ , which is much smaller than 
(Pc{d) ^ 2~'^. Therefore, the number of contacts between 
neighbors for p — Pc{d) could be large enough in order 
to have a cage. Furthermore there is evidence that the 
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packing fraction (pMRj{d) for maximally random jammed 
states in d dimensions is given by (ci + C2d)2~'^ [28|] or 
even with an additional term with a quadratic prefac- 
tor c^cP [29j. These densities are not larger than fc{d). 
The corresponding pair correlation function g2{T] d) flat- 
tens under an increase of d from 3 to 6 i.e. comes 
closer to the ideal gas value gig{r; d) = 1, for r > a. Sec- 
ond, concerning MCT S{k;d,(pc{d)) = 1 does not imply 
that the direct correlation function c{k; d, (pc{dj) and the 
vertices are zero. It is the quadratic dependence of the 
vertices on c(fc; d, tp) in combination with its explicit n- 
dependence and the use of the scaled variables k ^ ka/d 
which make the coupling of the modes finite, despite the 
small packing fraction (pc{d) ~ d^2~'^. This is completely 
similar to the MCT approach to colloidal gelation for 
a li quid of hard spheres with attractive Yukawa poten- 
tial [sol . IsTj . For packing fraction c/j — > and potential 
strength K ^ oo with Kip'^ — T — const, these authors 
prove that S{k) — > 1 for all k. However, the vertices re- 
main finite. At a critical value Tc a liquid-gel transition 
occurs. fc{k) is similar to fc'^\k), quite analogous to 
our outcome. The equilibrium structure at Tc is highly 
ramified where the "caging" of a sphere is generated by 
a smaller number of neighbors. 



C. Conclusions 

In section A. we have presented our various results 
from MCT in high dimensions. Now we want to discuss 
the most essential findings in the light of earlier results 
and will draw some conclusions. 

Our critical nonergodicity parameters have a non- 
Gaussian fc-dependenco, in variance with the assumption 
in Ref . [l^ . This discrepancy is the origin of the differ- 
ent pre-exponential factor of the critical packing fraction 
which we have found to be quadratic in d, and not linear 
[3l ■ Due to this quadratic d-dependence our MCT result 
for fcid) is larger than the Kauzmann packing fraction 
(fK (d) (Eq. ©). This cannot be true, since the packing 
fraction for the Kauzmann transition (static glass transi- 
tion) should be above the packing fraction for the MCT 
transition (dynamical glass transition). ipK (d) has been 
calculated within a small cage expansion (l8|, which is 
a kind of Gaussian approximation. This could be the 
reason why <fK (d) (Eq. (HD) is below ipc (d) (Eq. ((55)) ). 

As argued in subsection B. it is not necessarily true 
that a structureless static correlator rules out caging. 
But, even if the cage effect would be absent, the essential 
question would be whether the quality of both MCT ap- 
proximations necessarily requires the existence of caging, 
or not. Since an analytical investigation of the validity of 
these approximations seems to be extremely difficult, a 
way to get an insight is an approach by a computer simu- 
lation. Provided such simulational results would deviate 
more and more from our results with increasing dimen- 
sions this would hint at a failure of MCT for d — )■ oo. 
Such a failure would imply that MCT, which has been 



interpreted as a mean field theory does not become 
exact in the limit of high dimensions, in contrast to equi- 
librium phase transitions. We hope that these concluding 
remarks may stimulate and encourage further investiga- 
tions, contributing to a better understanding of MCT. 
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APPENDIX 

We want to prove the following statement: for all pack- 
ing fractions (p (d) such that 2'*(p (d) does not increase ex- 
ponentially or faster with d for d — > oo, the static three- 
point correlation function S^^'^{k, k') reduces in the limit 
of high dimensions to S (fc) S [k') S{\k + k'\), i.e. the con- 
volution approximation becomes exact. For this proof 
we use the Ornstein-Zernike equation for three-particle 
correlation functions 13211 : 



^(3) {k^k''^ = S {k) S [k') S{\k + k'\) (l + n^c^'^\k,k') 

. J48) 

So, we have to show that n^c''^^{k, k') for all fc, k' for 
d ^ oo and (p constrained as above. The explicit depen- 
dence of c^^^ {k, k') on k, k' does not have to be considered, 
as we can use for all k, k': 



,(3) 



{k, fc') 



d^.ydVe-e--c(3)(,.') 
< / d'^r [ d'^r' c^^) (f, r ') 



(49) 



Now we can expand c*-'^-' {r,r') into diagrams, where the 
lines are Mayer functions and the vertices are single par- 
ticle densities. This expansion only consists of loop di- 
agrams. We want to show now, that the contribution 



of each of these diagrams to J d'^r J d'' 



.(3) 



(f,r") 



vanishes in the limit d — )■ oo. To do so, we apply the theo- 
rem of Wyler, Rivier and Frisch 0. This theorem states, 
that a loop diagram leads to an exponentially smaller 
contribution to an integral like the one appearing in the 
last line of Eq. (|49p . than a tree diagram of the same or- 
der. The simplest diagram in the expansion of c*^"^) (r, r ') 
reads 



C(P (r, r') = (cr — r) G (cr — r') O (cr — |r — f' 



(50) 
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which can be inserted into Eq. (j49 



,(3) 



< 



(51) 



< d'^r d'^r' e (ct - r) e (cr - /) 9 (cr - |?^- r'\ 



The integral occuring in Eq. (|5ip leads to an exponen- 
tially smaller contribution than the corresponding tree 
diagram of the same order [2|]: 

d'^r j d'^r' 6 (cr - r) e (cr - r') 6 (cr - \r ~ r'\) < 

< a'^ j d'^r J d'^r' 6 (a - r) 9 (a - r') 

(52) 



where Vd (cr) is the volume of a d-dimensional sphere with 
radius cr and 



a < 1. 

From Eqs. ([?T|) and ([5^ we obtain: 

4'' (fc, k') I < a"^ (nT/d {a)y 



Together with 



^ = "^^d (I) = 2-''nFd (a) 



(53) 



(54) 



(55) 



this leads to 



(fc,fc')|<a'^ (2'^^)', 



(56) 



i.e. for all packing fractions, where 2'^(p does not increas- 
ing exponentially or faster with d, we obtain from Eqs. 



,(3) 



d^oo 



forallfc,fe'. 



(57) 



The contribution of all other diagrams are also exponen- 
tially smaller than the corresponding tree diagrams Q], 
which leads to: 



,(3) 



<afn^Vd{a)rn^ (VAa)Y (58) 



where j is the number of vertices over which it has to be 
integrated in the corresponding diagram and is always 
smaller than one. From this we obtain 



^^'^Ik.k' 



< Qfi (2 ip) 



(59) 



or 



for all fc, k' and all i (60) 



provided 2'^(p{d) does not increase exponentially or faster 
with d. 
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